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Abstract 

We show that colliding vortex beams instead of (approximate) plane waves can lead 
to a direct measurement of how the overall phase of the plane wave scattering amplitude 
changes with the scattering angle. Since vortex beams are coherent superpositions of 
plane waves with different momenta, their scattering amplitude receives contributions 
from plane wave amplitudes with distinct kinematics. These contributions interfere, 
leading to the measurement of their phase difference. Although interference exists for 
any generic wave packet collision, we show that using vortex beams dramatically enhances 
sensitivity to the phase in comparison with non-vortex beams. Since the overall phase is 
inaccessible in a plane wave collision, this measurement would be of great importance for 
a number of topics in hadronic physics, for example, meson production in the resonance 
region, physics of nucleon resonances, and small angle elastic hadron scattering. 

1 Introduction 

1.1 The phase of the scattering amplitude 

Scattering of energetic particles on targets and on each other is the standard way to inves- 
tigate their structure and interactions. Strength of a scattering process is characterized by 
the complex scattering amplitude M = \M\e ia . If we have in mind the usual situation of 
plane- wave scattering (initial and final state momenta are generically denoted as hi and kf, 
respectively), then it is only the absolute value of the amplitude \Ai(kf,ki)\ which is experi- 
mentally accessible via the cross section measurement, but not the overall phase a(kf, fcj). 

Despite not being directly measurable in plane wave scattering, a(kf,ki) is still a well- 
defined quantity which can be traced back to partial wave phase shifts. This phase can be 
of much importance, especially in hadronic physics, where a large phase is usually a sign of 
complicated dynamics. For example, in hadronic production or scattering experiments in the 
resonance region the phase of the amplitude traces a characteristic curve as a function of the 
invariant mass when one sweeps across a resonance. Another example is the high-energy elastic 
scattering of hadrons (notably, pp and pp). Due to the dynamical properties of the t-channel 



exchange with the vacuum quantum numbers, dubbed the Pomeron, (the hadronic part of) 
the scattering amplitude becomes mostly imaginary, but it also contains a real part. Their 
ratio, usually quantified by the parameter p(s,t) = ReA^/ImA^, turns out to be a sensitive 
probe of the Pomeron models, see for example [1] and references therein. 

Being so important for hadronic physics, the phase of the scattering amplitude is impossible 
or notoriously difficult to measure in the usual (approximate) plane wave scattering. Only 
when the reaction can proceed via the hadronic and the electromagnetic mechanisms and when 
the electromagnetic phase is sufficiently constrained, can the phase of the hadronic amplitude 
be accessed via the interference pattern between the two amplitudes. However, even this 
kind of measurement often poses experimental challenges. For example, in high-energy elastic 
processes the phase of the hadronic amplitude can be measured via the Coulomb-hadronic 
interference in a very small-i domain, [2], where the phase of the Coulomb amplitude itself is 
a matter of debates, see [3]. As a result, the p-parameter remains poorly constrained at high 
energies even for the pp/pp scattering, let alone other hadrons. The phase of the Pomeron 
amplitude in diffractive photo- or electroproduction of 7r + 7r~ pairs can be measured via its 
interference with the Primakoff mechanism, [I], but this suggestion has not yet been tested 
exp er iment ally. 

In other cases, when there is no clear reference mechanism with a known phase available, 
such as process 7p — > K + A, the overall phase determination becomes impossible, This rep- 
resents the notorious overall phase ambiguity in the interpretation of meson photoproduction 
data and the spectroscopy of baryonic resonances. Indeed, disentangling several interfering 
partial waves (and therefore the possible baryonic resonances in these channels) depends on 
how the phase changes with the polar angle of the produced meson. 

1.2 The essence of the proposal 

Here we propose a novel type of experiment which can probe the phase of the plane wave 
scattering amplitude a(kf,ki). To be more accurate, it measures how this phase changes 
with the production angle (its absolute magnitude still remains unmeasurable). This phase 
sensitivity comes from the fact that the scattering amplitude of two collinear Bessel vortex 
beams (to be defined below) is a coherent sum of two plane-wave amplitudes corresponding 
to fixed final momenta kf but slightly different initial momenta, ki and k[. These plane- wave 
amplitudes interfere, and as the result one gains sensitivity to their relative phase a(kf, ki) — 

In principle, interference among several plane wave amplitudes arises in scattering of any 
generic wave packets. The crucial feature of our suggestion is to exploit a very special form of 
transverse wave packets: the so-called vortex beams, in which wave fronts have helical shape 
and contain phase vortices. We will show that the phase sensitivity of the vortex beam collision 
is dramatically better than for any generic (for example, simple Gaussian) non-vortex wave 
packets. 

In order to make the proposed idea even more transparent, let us compare it with the 
method of partial waves, in which the initial and final states are expanded in spherical waves 
with the given orbital angular momentum I and its projection m. The phase at of the partial 
wave scattering amplitude is defined by its (real) phase shift 5e and the inelasticity of the 
scattering. Since the plane wave is a superposition of many partial waves, measurement of the 
full 4n angular dependence of the plane wave cross section can reveal the phase differences 
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ae — ai£>, but cannot access the overall plane wave phase a(kf,ki). On the other hand, if we 
were able to prepare and collide initial states as pure partial waves with given £, which would 
be superpositions of plane waves coming from all directions, we could measure individual 
cross sections og and deduce the angular dependence of the plane wave phase a(kf, fcj). The 
problem, of course, is that preparation of initial states coherent over the entire solid angle is 
unfeasible. 

A vortex beam is a cylindrical wave, which represents an intermediate case between a 
plane wave and a spherical harmonic. It is a superposition of plane waves with different but 
close momenta, and therefore vortex beam scattering is also sensitive to the relative phases of 
plane wave amplitudes. In contrast to spherical waves, however, vortex beams can be created 
experimentally. 

1.3 Vortex beams 

Optical vortex beams, that is, laser beams carrying non-zero orbital angular momentum 
(OAM) with respect to the beam axis, are well-known and routinely used in optics, [El Ej. 
Their application range from microscopy to astrophysics, |8j. Wavefronts of such a beam are 
not planes but helices, and each photon in this light field (a twisted photon) carries a well- 
defined OAM quantized in units of h. Vortex beams can exist for massive particles as well. 
Recently, following the suggestion of [9], electron beams carrying OAM were experimentally 
demonstrated, first using the phase plates [10] and then with the aid of diffraction gratings 
[TT] . Such electrons carried the energy as high as 300 keV and the orbital quantum number up 
tom~ 100. These vortex beams can be manipulated and focused just as the usual electron 
beams, and very recently remarkable focusing of a vortex electron beam to the focal spot of 
less than 1.2A in diameter was achieved |12j . 

It is conceivable that future progress will lead to creation of even more energetic twisted 
particles. For example, it was recently shown that X-ray beams carrying orbital angular 
momentum can be generated in ondulators and X-ray free-electron lasers, [13], and there 
exist plans to study them at the Argonne National Laboratory, [TJ]. Besides, it was noted in 
[T5] that Compton backscattering of twisted optical photons off an ultra-relativistic electron 
beam generates high-energy photons carrying non-zero OAM. Since the technique of Compton 
backscattering is well established pD] and is actively used in nuclear and particle physics (for 
example, at the SPring-8, [17], and HIgammaS, [15] . facilities), realization of this proposal 
seems feasible with today's technology. 

Vortex beams of other particles, for example of protons, can also be created experimentally. 
For low and intermediate energies one could use the same fork diffraction grating technique 
as used for electrons. Alternatively, one can exploit the vortex transfer which spontaneously 
happens in an elastic collision of a lower energy electron vortex beam with a high-energy 
proton beam that initially carried no phase vortex, [19]. 

Although these proposals have not yet been tested experimentally, one can start thinking 
about what new insights into properties of particles, and especially hadrons, one can gain in 
collisions of energetic vortex beams. First studies along these lines were conducted in [2D] , 
where it was noted there that such processes can probe observables which are not accessible 
in the plane wave collisions. The phase measurement discussed in this paper represents one 
particular example of such observables. 

The structure of the paper is the following. In Section [2] we introduce the formalism of 
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Bessel-type twisted states. Then, in Section [3] we develop further the theory of collision of two 
twisted particles and obtain the expression for its cross section. Section [4] explains how one 
can measure the angular dependence of the phase of the amplitude. In Section [5] we discuss 
what exactly is needed to actually realize the measurements we propose. Finally, in the last 
Section we draw our conclusions. 



2 Describing twisted states 

Here we briefly summarize the formalism of Bessel vortex states, which are the simplest states 
carrying a phase vortex. We focus on the scalar case and use the notation of [13] . We will 
accompany three-dimensional vectors with arrows, while the transverse vectors will be given 
in bold. 

A Bessel state is a non-plane wave solution of the free wave equation with a definite 
frequency u, longitudinal momentum k z , modulus of the transverse momentum |k| = x and 
a definite ^-projection of orbital angular momentum (orbital helicity) m. When written in 
cylindric coordinates r, <p r ,z, this solution |x, m) has the form 

pim<t> r 

| X; m ) = e*> t+tk * z ■ Vw(r) , V*m(r) = -=yfiiJ m (xr) , (1) 

v 2tt 



where J m (x) is the Bessel function. The transverse spatial distribution is normalized according 
to 

d 2 r ^*, m ,(r)Vw(r) = 5 mm > 8{>c - x) . (2) 

The azimuthal angle dependence oc e im< ^ r is the hallmark of the phase vortex. A twisted state 
can be represented as a superposition of plane waves: 

| X , m ) = e -**W J _0_ a _( k)e ^ j (3) 



where 



6(\k 



x 



axm (k) = (-i) m e %m *W2n Vl ' ' . (4) 



x 



This expansion can be inverted: 




i m e- imipk \>c,m) , x=|k|. (5) 



x 

m=— oo 



This invertible linear map between plane waves and twisted states means that twisted states 
simply represent another complete basis for transverse wave functions. When switching to 
vortex beams, we do not lose any asymptotic state we had before nor introduce any new state. 
More properties of twisted states can be found in [TSJ [20]. Here we only note that Eq. ^ in 
fact describes the passage from plane waves to twisted particles in description of a scattering 
process: one just needs to apply at the level of the scattering matrix the a^m- or a^ m -weighted 
transverse momentum integration for each incoming or outgoing twisted particle. 
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A pure Bessel state |x, m) with fixed x is non-normalizable in the transverse plane. A 
much better approximation to physically realizable states such as Bessel-Gaussian or aperture- 
limited beams is given by a fixed-m superposition of Bessel states 



with a properly normalized weight function /(x) peaked at x and having width a . This 
state, which we call the transverse Bessel wave packet, is normalizable (and localized) in the 
transverse plane. For calculational simplicity we assume that such states are monochromatic, 
that is, k z is supposed to vary with x so that k 2 z + x 2 is constant. Properties of such states 
were discussed in [21]. Here we note that they play the key role in resolving the non-forward- 
to-forward transition paradox in Bessel beam scattering which arose from contradicting results 
of [20] and PS]. 

When discussing twisted states of fields with polarization degrees of freedom, extra care 
should be taken. For example, in the case of photon the issues arising span from the absence of 
a well-defined polarization vector, which is replaced by the polarization field, to the problem 
of gauge-invariant separation of spin and OAM operators. Fortunately, these problems can be 
avoided in the paraxial approximation, namely at x < w, where the polarization and OAM 
degrees of freedom can be reliably separated, see discussion in [22]. The polarization field can 
be then approximated by a constant polarization vector lying in the plane orthogonal to the 
axis z. The case of twisted fermions, for which the paraxial approximation is also assumed, 
can be treated similarly. A detailed account of the ultrarelativistic twisted electron beyond 
the paraxial approximation was recently given in [23J. 

3 Beam-beam scattering beyond the plane wave ap- 
proximation 

3.1 General description 

The general theory of scattering of arbitrarily shaped beams was developed in [24J. Its for- 
mulation took into account not only the arbitrary spatial distribution, but also a non-uniform 
temporal dependence of the non-monochromatic beams. The basic objects describing the 
beams were not the coordinate or momentum wave functions, but the Wigner distributions 
which arose after statistical averaging over the ensemble of particles in each beam. This formal- 
ism was especially useful for description of coherent processes such as beamstrahlung, in which 
a particle from one beam coherently scatters from many particles in the counterpropagating 
beam. 

For the purpose of the present paper (namely, to demonstrate sensitivity of the vortex 
beam scattering to the phase of the amplitude) such a general set-up is not needed. Here, we 
make the following assumptions: (1) the colliding particles are in pure states describable by 
wave functions, (2) the beams are monochromatic. These two assumptions are made solely to 
simplify the calculations; if needed, one can redo the analysis in the most general framework of 
[23]. Finally, at a later stage we will also use the fact that the colliding beams are paraxial, and 
the longitudinal distribution of their envelopes is approximately uniform over a long distance 
compared to the transverse size. 




(6) 
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We consider a 2 — > 2 scattering (not necessarily elastic) and start with the plane-wave 
collision when initial particles with four-momenta k\ and hi scatter into final particles with 
momenta k[ and k' 2 . The total momentum is denoted as K = k± + k 2 = k[ + k' 2 . The ^-matrix 
element of this process is 

S PW = l {2nf^\k 1 + k 2 - k[ - K)^fyt§= , (7) 

where the invariant amplitude Ai is calculated by the standard Feynman rules. This ampli- 
tude, in general, depends on all the momenta involved, but for future convenience we explicitly 
indicate only the initial momenta. 

Let us now assume that the initial particles are not plane waves but are represented by 
coordinate wave functions ipi (f) and ip 2 if) normalized by J d 3 r\ifji(r)\ 2 = 1. The corresponding 
momentum-space wave functions are 

V(k) = J rf 3 r^(f)e-^, j ^- 3 \<P$)\ 2 = 1- ( 8 ) 

If these wave functions are not normalizable in the entire space, then a large but finite inte- 
gration volume is assumed in order to provide the regularization. 

The S"-matrix element for scattering of this initial state into the same plane-wave final 
state with momenta k[ and k' 2 , can be written as 

S = J J^^^ kl) ^ k2)Spw • (9) 

Since the beams are monochromatic, the number of scattering events into a given differential 
volume of the final phase space per unit time is 

(2nYS(E) d% d% 

\E X E 2 1 1 (2vr) 3 2^ (2tt) 3 2£ 2 ' 1 1 

Here 5(E) stands for 8(E X + E 2 -E[- E' 2 ), and 

ci 3 fci d 3 k 2 



(2tt) 3 (2tt) 3 



^i(fci)v?2(fc 2 )<5 (3) (^i + h - K)M(h, k 2 ) . (11) 



Let us now discuss the obvious but important issue of the final momentum distribution of 
the produced particles. In the conventional plane-wave case the final momenta are constrained 
by the overall momentum conservation. The value of K is then fixed, and the momenta k[ 
and k' 2 become maximally correlated. The usual procedure then is to use the momentum 



delta-function and remove, for example, the kl 2 integration. In our case, (10), this maximal 
correlation is absent. One must consider distribution of the events over momenta k[ and k' 2 , 
or alternatively, over k[ and K, where the exact K region depends on tpi(ki). 

3.2 Approaching the plane wave limit 



Before proceeding further with vortex beams, let us first see how (10) simplifies in the plane 



wave limit. The plane- wave limit corresponds to very compact momentum wave functions 
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(fii(ki) localized near (hi). The matrix element can then be approximated as Ai(k±,k 2 ) ~ 
A4((ki), (k 2 )) = Aio, and the expression for / becomes 

F = M j ^- 3 e-^Mr)Mr) • (12) 
Changing d 3 k' x d 3 k' 2 to d 3 k[d 3 K and integrating \ f\ 2 over K, one gets 

d 3 K\ff = \M \ 2 J £^\Mr)\ 2 \Mrl\ 2 ■ (13) 
This means that in the plane wave limit, ki = (ki), one can effectively replace 

|F| 2 \M,\^ 3 \h + k 2 -K)J ^y 3 \M^\ 2 \M^)\ 2 • (14) 

The number of events can therefore be split into the cross section and luminosity factors: 

du = da ■ L , (15) 

_ {2ixY5\k l + k 2 - k[ - k' 2 ) 2 d 3 k[ d 3 k' 2 
~ AE^v 1 ol {2^) 3 2E[ (2tc) 3 2E' 2 ' 

L = v J d 3 r n 1 (r)n 2 (r) , n,i{r) = \tpi(r)\ 2 . (16) 

Note that we inserted here by hand the relative velocity of the two plane waves, v = \v\ — v 2 \. 

It must be stressed that separation of the number of events into the differential cross 
section and the (conventional) luminosity is uniquely defined only for plane waves. Extending 
this splitting for non-plane-wave collisions is a matter of convention. In this way, one needs 
to (somewhat arbitrarily) introduce the notion of generalized cross section, [21], for example, 



by dividing the full du in (10) by L (16) with v defined for (ki) rather than ki. With this 



definition, the generalized cross section for non-plane-wave scattering takes form 



(2n) 3 \F\ 

da = da Rd 3 K, R = \ > 1 , (17) 

|yVIo| J d A r ni{r)n 2 {r) 

where do~o is the plane-wave .^-integrated cross section. In this notation, the plane-wave limit 
corresponds to R — > d^^ki + k 2 — K). 



3.3 Pure Bessel states 

Let us now consider the particular case of pure Bessel states for both initial particles. This 
problem was first considered in [20J, where this process was termed the "double-twisted scat- 
tering" and a generic expression for the cross section was found. That expression used a 
definition of the flux factor (i.e. luminosity) different from what was described above, there- 
fore that cross section differs in the overall normalization from our results. However, the main 
finding remains unchanged: for each K there exist two kinematic configurations which add up 
coherently. Since this finding is the key to our proposal, we discuss it below in more detail. 
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Consider a collision of two Bessel states |xi,mi) and |x2,m2}. For simplicity we assume 
that both vortex states are coaxial, that is, they are defined with respect to the same axis z. 
The wave functions describing these states are given by regularized versions of Q and (§. 
Namely, we assume a large but finite cylindrical quantization volume with macroscopic length 
L 7 and radius R. The coordinate wave functions are then 



y Zn \ 2n ~ 

where the orbital helicity m 2 is accompanied by the sign minus because the second particle 
moves in the negative z direction. The normalization coefficient c is fixed by 

1 = \c\*L M Xi [ rdrJ^iw) « |c| 2 ^, (19) 
Jo 71 

where we used the large-argument asymptotics for the Bessel function. The wave function in 
the momentum representation, for example, of the first particle is 

(pi(k) = 2nc5 r (k z - k u ) ■ (-i) m e imi ^J—5 r (Ikl - xi) , (20) 

V 

where two 5 r -functions denote regularized versions of the delta-function. For the longitudi- 
nal/transverse dynamics, 5 r .(q) is localized within |g| < 1/L Z and |g| < 1/R, respectively, and 
its integral over all q gives unity. The longitudinal factors can be easily accounted for, and 



the cross section (17) can now be represented as 



- da ■ R± , where R ± = — ..^i ^1 , , \ , (21) 



d 2 K u ^' ^ |M | 2 /d 2 rni(r)n 2 (r) 

where 

F± = J ^^a J , 1 , mi (k 1 )a^,_ m2 (k 2 )5( 2 )(k 1 + k 2 -K)>l. (22) 

Although we refer here to a Xjm (k) defined by Q, it is assumed that the delta-function inside 
it is the regularized one. 




Figure 1: Two interfering kinematical configurations on the transverse plane for vectors k x 
and k 2 of fixed length which sum up to the same vector K. 
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The integrals of this type in the limit R — » oo were studied in [20] ■ In this limit, a^m's 
contain true delta- functions, so that all four integrations in (22) are eliminated. It turns out 
that this integral is non-zero, only if there exists a triangle with sides x 1; x 2 , |K|, so that K 
lies in the ring 

|xi — x 2 | < |K| < xi + X2 . (23) 

For any K lying in this ring there exist exactly two kinematic configurations shown in Fig. [T] 
for ki and k 2 with azimuthal angles 

/x 2 2 + K 2 - x^X 

01 = <p K ± 61 , 02 = 0A' ± (-^2) , h,2 = arCCOS ( ' 2 ^ pj^j — - I . (24) 

Note that by definition < 5i + <5 2 < 7r. The momenta ki and k2 corresponding to the sign + 
in both ±'s are denoted as k+ and k^ (Fig. [TJ left), while for the sign minus they are denoted 
as k^~ and k 2 (Fig. [TJ right). The amplitude Ai taken at these two pairs of momenta will be 
denoted by Ai+ and Ai-, respectively. With this notation, the integral F± can be written as 

_ ( ^ m l~ m 2 c i(mi-m 2 )4>K V XlX2 I ^ e im 1 8 1 +im 2 S2j^^ _|_ e -«mi5i-im 2 <52^ ^ ^25) 

Here 

1 



A = _y 2xfx| + 2xfK 2 + 2x|K 2 - xf - x| - K 4 = -xix 2 sin(5! + 8 2 ) , (26) 
is the area of the triangle with side xi, x 2 , |K|. We then obtain 

2 XiX 2 



(2tt) 6 4A 2 



[\M + \ 2 + |7W_| 2 + 2Re (e 2i(mi5l+m2S2) M + M*_)} . (27) 



The last term here is the interference term (or the autocorrelation function), which is 
sensitive to the phase difference of the amplitudes M. + and M._. It is this term that allows 
one to measure how the overall phase of the amplitude changes with angles, which will be the 
subject of the next Section. 



3.4 End point singularities 



| 2 cannot be integrated, as it stands, in the entire 
K|, where A = 0, it develops singularities. Indeed, 



It is clear that the expression (27) for \F_ 
|K| domain. At the end points xi = X2 ± 
at these points M. + = Ai-, while 5i = or tt, so that the expression in brackets has a finite 
limit. Changing the integration variable from K 2 to a = 5\ + we get a non-integrable 
singularity of the type 



dK 2 



v , oc / — . (28) 
A^ ,/ sin a 

This singularity can be traced back to using the R — > oo limit. However, in this limit the 



da 



denominator of R± (21) also diverges, and one gets the indeterminate form oo/oo. The remedy 
in this situation is to consider large but finite R and take the limit of the ratio R± rather 
than its numerator and denominator separately. Then, the exact delta-functions inside a^m's 
in ((22l are replaced by regularized functions 5 r ,. The result (27) stays valid until |K| comes in 



the 1/R vicinity of the end points, where the function F± gets regularized. Therefore, at finite 



9 



R the effective integration limits in (28) are «o and tt — ao with ao ~ 0(1/ xR), where in the 
logarithmic accuracy x can be taken equal to any of Xj. After integration, the numerator of R± 
contains a large logarithm log(xi2). The denominator of R± (21) is also enhanced by a similar 



logarithm. Indeed, for large R the integral J d 2 rii(r)n 2 (r) is approximately proportional to 



R 



dr 



cos 



X\T 



7T7TT.1 



2 / irm 2 
cos ( x 2 r H — 



oc log(xi?) . 



Therefore, J d 2 kR± stays finite in the R — > oo limit. 

Following this analysis, one is tempted to conclude that in the R — > oo limit the only 
contribution comes from the two end points in the |K| region, while the contribution from 
the majority of the interval (23) is suppressed as l/log(xi?) and eventually vanishes. The 
tricky thing, however, is that this suppression is not so strong. Taking any reasonable macro- 
scopic R and microscopic x, one gets log(xi?) ~ 50. So, the inner region of the |K|-interval 
does give a non-vanishing contribution even for macroscopic quantization volumes. This is a 
curious example of a situation when the unphysical limit of infinite quantization volume is 
quantitatively different from the physically motivated extremely large but finite volume. 

The case of a pure Bessel state, although important for the formalism itself, is of academic 
interest for applications. Indeed, true Bessel beams of energetic particles with a macroscopic 
spatial extent cannot be created experimentally. Taking into account the vortex beams already 
created and envisaging the suggestions for generation of high-energy particles with OAM 
[T5"j [T9"] . one can count only on transverse Bessel wave packets ^ with the weight function 
/(x) having a reasonably narrow but finite width a . In this case, the ratio R± contains the 
square of 

F\ — [ dx 1 dx 2 f 1 (x 1 )f 2 (x 2 )F ± (29) 



instead of just F±. In the case a <C x, one can basically repeat the above analysis using the 
regularization parameter 1/a instead of R. In this case, the logarithmic factor log(x/<r) is not 
big at all, and the |K| integration gives comparable contributions from the end points and 
from the central part of the |K|-region. 



0.8 



0.4 





m, = 0, m 2 = 




\ a i = K i 1 50 / 






a. _ k j / 5 




: 











- 



0.5 



0.4 



0.3 - 



0.2 



0.1 - 







= 3, m 2 = 










Oj = k„ , / 50 










a . - k j / 5 






u 








l 



0.4 



0.6 



1.2 



1.4 



1.6 
K 



0.4 



0.6 O.f 



1.2 1.4 



1.6 
K 



Figure 2: Quantity |K|i?j_ as a function of |K| (measured in arbitrary units) for mi = (left 
pane) and m\ = 3 (right pane); m 2 = in both cases. The plots were drawn for xi = 1 
and x 2 = 0.5 measured in the same arbitrary units. The solid and dashed lines correspond to 
&i = ^oi/50 and <7j = x 0i /5, respectively. 
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For illustration, we show in Fig. [2] the ratio R± multiplied by |K| for the simplest case 
of constant scattering amplitude: Ai + = M._ = Ai . The qualitative dependences discussed 
above are well reproduced. For sufficiently small (7j this quantity indeed peaks near the end 
points (which are at 0.5 and 1.5 in the units chosen), so that the end points indeed dominate, 
especially at non-zero orbital helicities. For not too small a, such as <7j = xoi/5, the end points 
do not dominate over the inner region at all; however, in this case the non-zero helicities become 
noticeably suppressed. 



4 Measuring the phase of the amplitude 

It was already mentioned above that the expression for the (generalized) cross section in a 
vortex-beam collision is sensitive to the difference of the phases of the invariant amplitudes 
calculated at two sets of the initial momenta. This sensitivity comes from the interference 



term in (27). In this Section we quantify this sensitivity. 

Although the results derived below are rather universal and not particle-type specific, we 
will use for convenience a notation that alludes to one particularly representative case: the pion 
photo-production 7^ — > up. In this situation, we will refer to the particles with momenta k\ 
and ki as the photon and proton, respectively, while k[ will denote the final pion momentum. 

Since the purpose of this paper is to demonstrate the sensitivity to the phase in a generic 
set-up, we will not stick to any particular model for this specific process. Instead we consider 
a toy model, in which the modulus of the amplitude is constant while the phase depends only 
on the angle between the photon and the pion in the proton rest frame, 6* n : 

M = \M\e ia ^. (30) 

The assumption that is constant is not, in fact, important, as it will be discussed at the 
end of this Section. 



4.1 Differential cross section 




Figure 3: Schematic representation of the pion photoproduction kinematics in the proton "rest 
frame" in the case when both the photon and the initial proton are represented by the Bessel 
vortex states |x 7 ,m 7 ) and |x p ,m p ). 

Our goal now is to find an efficient way to extract phase a (or, to be more precise, a' = 
da/d9*) from the vortex-beam photon-proton cross section. To this end, we consider the 
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double-twisted 7p collision, in which both photon and proton are represented by Bessel beams 
(which will later be replaced by the transverse Bessel wave packets) with orbital helicities 
m 7 and m p . The two twisted states are defined with respect to the same z axis, so that a 
longitudinal boost along this axis does not change the (transverse) vortex structure of these 
states. We are allowed, therefore, we switch to the frame where the longitudinal momentum 
of the proton is zero: this is the closest approximation to the rest frame one can achieve with 
twisted particles. The kinematics of this process is shown in Fig. |3j 

The differential cross section depends on \F±\ 2 given by (27), and in our toy model it is 
proportional to 

— oc 1 + cos [2(m 7 5 7 + m p 5 p ) + (a + - a_)] , (31) 

where the phases a± correspond to the two kinematical configurations shown in Fig. [TJ 

Suppose that the spherical angles of the initial photon are 8 1 <C 1, 7 , and the angles of 
the final pion are 8 n , <p^. In order to calculate the phases a±, one has to take into account 
two kinematical effects. First, the angle between the photon and the pion differs from 6 n : 







77T 



T , 



# 7 cos(0 7r — 



(32) 



Second, in the frame chosen the collision involves a proton which is not at rest but moves with 
a transverse momentum k p , |k p | = x p . Since the angle 8* n is defined in the proton rest frame, 
one needs to perform a transverse boost with speed v = x p /E p , where E p w M p . Keeping 
only terms which are linear in this velocity, we can relate the angles 6* and 6 yn with the 
pion energy and momentum in the lab frame (E n , p n ) and in the proton rest frame (E*, p*): 
E n - p w cos # 77r = E* — p* cos # 77r . Then, we get 



8* w ,/ 



8~fTT ~{ 



Mr, 



1 - 



COS^ 



COS (0TT — 



(33) 



Together with (|32|), this equation finally gives 8*^. 

f»i and <f) p = 02 for the two kinematical configurations are given 



Now, the values of 7 



by (|24|). The difference between the two phases is then 

E~, 



a_i_ — a_ 



a 



-2a' 8, 



1 + 



M„ 



cos8 n 

~J7 



sin 5 y sin(<p n — <fi K ) , (34) 



where M p is the proton mass and we used here x 7 sin<5 7 = x p sm5 p . Thus, we finally obtain 
the expression for the K-differential cross section: 



da 
d?K 



oc 1 + cos[2(m 7 <5 7 + m p 5 p )] + 2a' 8^ 



1 + 



Ej 



COS 8-n 



x sin[2(m 7 <5 7 + m p S p )} sin <5 7 sin(0^ — 4>k) ■ 



(35) 



Note that this expression depends both on modulus |K| via 8 1 , 5 P , and on the azimuthal angle 
6 K . 



12 



4.2 Azimuthal asymmetry 



Our goal now is to extract the last term in (35) which contains the desired quantity a'. Before 



we proceed, let us note the importance of using vortex beams and not just arbitrary wave 
packets. Indeed, if both initial particles did not have any phase vortex, then m 7 = m p = 0, 
and the last term simply vanishes. In order to get the sensitivity to the phase, one would 
need to consider the next term in the Taylor series, which would be quadratically suppressed 
at small # 7 . It is the non-zero orbital helicity that allows us to retain the term linear in small 
# 7 and hugely boost the sensitivity to the phase. 



Next, note that if we integrated (35) over all (fix, the last term would vanish. In order to 



extract it, we introduce the sine- weighted cross section 

/d 

to which only the last term contributes. Of course, this definition is not unique, and one 
could as well use the sign weight factor sign[sin(0 7r — 4>k)] instead of the sine itself. The key 
observable is then the azimuthal asymmetry: 

A=—, where a = [ d 2 K^y . (37) 
a J crK 

To avoid misunderstanding, we remind the reader that both quantities, Act and a, are still 
differential in dQ n . 




Figure 4: Detecting non-collinearity of the two final transverse momenta is a prerequisite to 
measuring the azimuthal asymmetry. 



Let us stress that in order to be able to extract the azimuthal asymmetry (37) from the 



experiment, one must be able to measure the transverse momenta of the pion and the final 
proton with a sufficient precision so that the non-collinearity of k^ and k p / can be observed. 
Basically, one must be able to detect the sign of sin(0 7r — 4>k), which is illustrated in Fig. |ij 
Any non-zero asymmetry measured in the experiment will be an indication of a non- 



constant phase. Using (35), one can write the azimuthal asymmetry as 



A = a' ■ P , 

where the analyzing power P can be written as 

M p V At 



V 



(38) 



(39) 
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The quantity V for pure Bessel beams can be represented as the ratio 



V 



dK 2 s ^ n (^ m 7^'r + 2m p 5 p ) sin<5 7 



sin (£ 7 + 8 P ) 



2 1 + cos(2m 7 5 7 + 2m p 5 p ) 



(40) 



sin (<5 7 + 5p) 



If instead of pure Bessel beams one uses transverse Bessel wave packets, then the general 
expressions for the azimuthal asymmetry (37) and the analyzing power (39) remain the same 
(with replacement 6* 7 — > 9 0l = x 07 /i? 7 ), but the expression for the quantity V becomes more 
complicated: 



( dK 2 F 1 F 2 
V = r „ rn „ n , where 



2 ' 



/ dK 2 F{ 



Fi= dK^dH p f J (^)fp(x p 



dx^dxpf^x^jfp^p 



(41) 



cos(m 7 <5 7 + m p 5p) 
A /x 7 x p sin(5 7 + 5 P ) 

K~ t sin 5-y sin(m 7 5 7 + m B 5, 



p w pj 



xq 7 v /x 7 x p sin(5 7 + 5 P ) 



(42) 




Figure 5: The ratio V defined in (41) for the transverse Bessel wave packets as a function of 
photon's orbital helicity m 7 ; the proton's orbital helicity is fixed at m p = 0. The calculations 
are done for x 7 = 2x p ; the solid and dashed lines correspond to Oi = x i/5 and Oi = x j/50, 
respectively. 



Let us discuss some properties of the analyzing power P, (39). First, it is linearly sup 



pressed by the small opening angle of the photon Bessel state, 9 1 = x 7 /_E' 7 1. The factor 
in brackets is of order one; boosting it further by increasing the photon energy does not help 
much because 9 1 will simultaneously decrease. Finally, the last factor, V, shows basically the 
relative strength of the sine and cosine of the angle m 1 5 1 + mp8 p averaged with the appropriate 
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weight functions. Of course, the denominator, as it stands in (40), is divergent due to the end 



point singularities discussed in Section |3.4| On the other hand, the numerator is convergent 
due to the extra sin<5 7 factor which turns zero at the end points. Thus, for the pure Bessel 
states the quantity V is logarithmically suppressed. However, for the transverse Bessel wave 
packets this suppression is relieved. In Fig. [5] we show this quantity for zero proton orbital 
helicity for a range of photon orbital helicities m 7 . One sees that for not too large m 7 this 
quantity is of order one. In can be checked in addition that is if x 7 <C x p or x 7 3> x p , the 
quantity V is additionally suppressed; therefore, x 7 ~ x p will be the optimal choice. 

Let us now summarize this analysis in a list of rules which one should follow in order to 
maximize the analyzing power. When preparing the vortex states of the photon and proton, 
one should (1) select small but non-zero orbital helicities, (2) choose x 7 ~ x p , (3) choose not 
too small o"j/xoi. If these criteria are met, one can expect the azimuthal asymmetry of the 
order 

A ~ a' 9 7 . (43) 

This estimate can be safely used in qualitative discussions of the effect. 

Let us finally discuss whether the assumption of the constant absolute value of the ampli- 
tude was essential in our analysis. In fact, it was not. Indeed, even if \A4\ changes with 
angles but phase a remains constant, then expression in the brackets in (22) can be written 

as 

e ia [cos(m 7 5 7 + m p 5 p )(\M + \ + \M-\) + isin(m 7 5 7 + m p 5 p )(\M+\ - \M-\)] . (44) 

The antisymmetric factor sin(0 7r — </>#-) appears only in the difference of the absolute values. 
Since there is no interference between the two terms here, \F±\ 2 will contain sin 2 (0 7r — 4>k), 
and the azimuthal asymmetry will be absent. Considering transverse wave packets does not 
break the relative 7r/2 phase between these two terms and does not change this conclusion. 
Therefore, even if is not constant, one can safely associate non-zero azimuthal asymmetry 
with varying phase. 



5 How far is this suggestion from experimental realiza- 
tion? 

The measurements suggested in this paper exploit the non-plane-wave nature of the colliding 
beams. Although it is not surprising that a wave packet collision is kinematically different 
from the plane wave collision, we would like to stress once again the unique advantage of using 
vortex beams. If both colliding beams were regular wave packets without phase vortices, then 
the phase of the amplitude would be effectively averaged over a certain interval of 8* /7T with 
a small width # 7 . The difference between this average (a) and its plane-wave value would be 
quadratically suppressed, 0(9T). In our case, using vortex beams with a reasonably narrow 
x-distribution and measuring the azimuthal asymmetry, one simply removes the central value 
of a and has a very clean phase difference between two kinematical configurations. This 
difference is only linearly suppressed by the small 9 T 

Next, it is true that the experiments proposed here are not yet feasible with today's 
technology. As we mentioned in the introduction, so far only twisted electrons and photons 
with low to intermediate energies and small values of x have been created experimentally. 
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Though there exist suggestions of how to bring vortex beams into high-energy physics, their 
realization requires dedicated accelerator and beam physics studies. 

To be more specific, the following three objectives must be achieved in order to make our 
proposal feasible. 

• Both colliding particles must be wave packets of Bessel states, and at least one of them 
must contain a phase vortex, that is, must carry orbital angular momentum. This 
is necessary if one wants an observable which is only linearly and not quadratically 
suppressed by small B 1 . Arguably, the simplest way would be to create energetic twisted 
photons, while for the protons it is sufficient to take a tightly focused Gaussian wave 
packet without any phase vortex, which corresponds to the transverse Bessel wave packet 
with m p = and a v ~ x p . 

• The transverse momenta inside the Bessel states, Xj, must be sufficiently large, namely, 
larger than the transverse momentum resolution of the detector. This requirement makes 
it possible to detect non-collinearity of the transverse momenta of the final particles, to 
determine the value of sin(0 7r — <p K ) or at least its sign and, therefore, to reconstruct 
the azimuthal asymmetry. Besides, large values of x 7 are needed to keep the analyzing 
power not too small. 

• The initial vortex state must be prepared with non-zero but small orbital helicity, with 
x 7 ~ x p and not too small Cj/xj. 

Among these requirements, arguable the most difficult to satisfy is the second one. For a 
typical hadronic experiment, it implies x ~ O(10 — 100 MeV), which corresponds to focusing 
an incoming particle to a spot of about 10 fm. What has been achieved experimentally so far 
is x ~ 10 keV and focussing to lA. Increasing x by at least three orders of magnitude seems 
to be a real challenge. Nevertheless, the huge expected payoff makes it worth trying to push 
towards these objectives. 



It is important to note that theoretical predictions of the azimuthal asymmetry (37) depend 
not only on the underlying microscopic physics, but also on the exact shape of the colliding 
vortex beams. However detection of any non-zero azimuthal asymmetry will signal the phase 
variation with the angles, as it cannot be mimicked by the absolute value of the amplitude. 



6 Conclusions 

In this paper we argued one can measure the angular or momentum transfer dependence of 
the overall phase of the plane wave scattering amplitude by colliding vortex beams. This 
measurement would be of much importance in several subfields within hadronic physics. 

Since the physics of energetic vortex beams is only emerging, the proposed experiments 
cannot yet be realized with today's technology, and further research in the beam physics is 
needed. We discussed the requirements that must be met in order for the measurements 
proposed to be feasible. In general, we see this proposal as an example of how beneficial 
vortex beam collisions can be to high-energy physics as they can probe quantities which are 
inaccessible in plane-wave collisions. 
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